Gravity and Large-Scale Non-local Bias 



Kwan Chuen Chan, Roman Scoccimarro, 1 and Ravi K. Sheth 2 ' 3 

1 Center for Cosmology and Particle Physics, Department of Physics, 
New York University, NY 10003, New York, USA 
2 Center for Particle Cosmology, Department of Physics and Astronomy, 
University of Pennsylvania, 209 S 33rd Street, Philadelphia, PA 19104, USA 
3 Abdus Salam International Centre for Theoretical Physics, Strada Costiera 11, 34-151, Trieste, Italy 

(Dated: February 22, 2012) 

For Gaussian primordial fluctuations the relationship between galaxy and matter overdensities, 
bias, is most often assumed to be local at the time of observation in the large-scale limit. This 
hypothesis is however unstable under time evolution, we provide proofs under several (increasingly 
more realistic) sets of assumptions. In the simplest toy model galaxies are created locally and linearly 
biased at a single formation time, and subsequently move with the dark matter (no velocity bias) 
conserving their comoving number density (no merging). We show that, after this formation time, 
the bias becomes unavoidably non-local and non-linear at large scales. We identify the non-local 
gravitationally induced fields in which the galaxy overdensity can be expanded, showing that they 
can be constructed out of the invariants of the deformation tensor (Galileons) , the main signature 
of which is a quadrupole field in second-order perturbation theory. In addition, we show that this 
result persists if we include an arbitrary evolution of the comoving number density of tracers. We 
then include velocity bias, and show that new contributions appear; these are related to the breaking 
of Galilean invariance of the bias relation, a dipole field being the signature at second order. 

We test these predictions by studying the dependence of halo overdensities in cells of fixed dark 
matter density: measurements in simulations show that departures from the mean bias relation 
are strongly correlated with the non-local gravitationally induced fields identified by our formalism, 
suggesting that the halo distribution at the present time is indeed more closely related to the mass 
distribution at an earlier rather than present time. However, the non-locality seen in the simulations 
is not fully captured by assuming local bias in Lagrangian space. The effects on non-local bias seen 
in the simulations are most important for the most biased halos, as expected from our predictions. 

Accounting for these effects when modeling galaxy bias is essential for correctly describing the 
dependence on triangle shape of the galaxy bispectrum, and hence constraining cosmological param- 
eters and primordial non-Gaussianity. We show that using our formalism we remove an important 
systematic in the determination of bias parameters from the galaxy bispectrum, particularly for 
luminous galaxies. 



I. INTRODUCTION 

Galaxies are one of the main probes of modern cosmol- 
ogy. However, the galaxy clustering amplitude depends 
on galaxy type, so not all types can be unbiased trac- 
ers of the dark matter pQ. Therefore, understanding and 
accounting for this bias is important. It is common to as- 
sume that this bias is a local and deterministic function 
of the dark matter density field (e.g. other properties of 
the field than the local overdensity, such as the tidal field, 
are assumed to produce negligible effects on the galaxy 
distribution) , so the galaxy density contrast at any given 
time can be written as a Taylor series in the dark matter 
density at that time (2J. One of the main goals of this 
paper is to show that, if there is any time at which this 
is a good approximation, then it is not good at any other 
time. A related goal is to argue is that this should be a 
better model at early than at later times, in a sense that 
will be made more precise later in this paper. 

In the galaxy distribution, one expects departures from 
the local deterministic bias model on scales where non- 



linear baryon physics matters. Nonetheless, on the scales 
larger than those associated with galaxy formation pro- 
cesses, the deterministic local bias is expected to be ac- 
curate, except for a possible constant shot-noise-type 
contribution [3J H]. This has motivated the use of the 
deterministic local biasing prescription for interpreting 
clustering measurements in galaxy surveys, in particular 
this model has been heavily used in interpreting mea- 
surements of three-point functions and other measures of 
non-Gaussianity [2] l5l-f2T]. However, it is also common 
to assume that galaxies are closely associated with dark 
matter halos [22]. So it is natural to ask if halo bias 
is a deterministic function of the local dark matter over- 
density. Numerical simulations indicate that, on scales of 
order 20 Mpc and less, halo abundance is not a determin- 
istic function of the dark mass [251125] . This manifests 
as stochasticity in the relation between the galaxy and 
dark matter density fields at the present time (26] [27] . If 
we distinguish between the stochasticity associated with 
some initial or formation time, and that due to evolution 
from this time to the time of observation - then the ques- 
tion arises as to which matters more on the large scales 
which the next generation of galaxy surveys will probe. 
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In the Excursion set model of halo formation, abun- 
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dance and clustering, it is the initial fluctuation field 
which is fundamental [28- 30 . In these models, the origin 
of the first source of stochasticity is relatively straight- 
forward to understand: the initial random fluctuation 
field is expected to have structure on arbitrary small 
scales, so the substructure within large patches of the 
same large scale overdensity may differ from one patch 
to another. Whether or not a small patch forms a halo 
is known to be closely related to the initial overdensity 
of the patch. If the initial overdensity is the only param- 
eter which matters (e.g., in spherical evolution models, 
where tidal fields etc. play no role) then the fact that 
the small scale density is correlated with the density on 
larger scales produces stochasticity in halo abundances 
within large spheres of fixed initial overdensity. Much of 
this scatter is just a sort of shot-noise which decreases 
as the cell size is increased [24 . So, on large scales, a 
deterministic model for the bias can be quite accurate. 
If halo formation depends on quantities other than local 
density [3H 131] , then this may contribute to the stochas- 
ticity seen in the initial conditions. But if these other 
quantities are correlated over shorter scales than is the 
density, then their effects will be subdominant on large 
scales, and so they may be neglected in studies of suffi- 
ciently large scale bias. In what follows, we will assume 
this is the case. 

That is to say, the main goal of this paper is to study 
departures from the local deterministic bias model which 
may appear on scales larger than those associated with 
galaxy or halo formation (i.e., above a few tens of Mpc). 
We will show that, even if the bias is local and deter- 
ministic at some given time (which we will usually call 
the formation time), then subsequent nonlinear gravita- 
tional evolution will generate non-local bias. In this re- 
spect, our results serve as a well-motivated model for 
non-local bias. Other works on non-local bias have pro- 
vided models [H 1331 131] based on statistical (as opposed 
to dynamical) considerations. The virtue of our approach 
is that it gives a concrete form of non-local bias that must 
be present even if formation bias is truly local, and we 
demonstrate for the first time their presence in numerical 
simulations. In addition, we show that our non-local bias 
model solves a systematic effect in the determination of 
the linear bias from bispectrum measurements for biased 
tracers. 

Since evolution plays an important role in the discus- 
sion, we devote a substantial part of this paper to the 
study of the evolution of bias and how it generates non- 
local bias. The evolution of halo bias, under the assump- 
tion that the number of halos was conserved and their 
motions were not biased relative to the mass, was first 
studied by [23]. They showed that the predictions for 
this evolution, based on a spherical collapse model for the 
dynamics, provided a good description of how halo bias 
evolves. At linear order (linear theory evolution of the 
linear bias factor) , this calculation agrees with that from 
combining the continuity equation with perturbation the- 
ory, again assuming no velocity bias [35] . At linear order, 



the perturbation theory approach can be generalized to 
include stochasticity and galaxy formation as a source 
[361 137] . However, going beyond linear order, either in 
evolution or in bias, is less straightforward. 

Evolution of the higher order bias factors was investi- 
gated in [33J [3S] , but these works approximated the non- 
linear gravitational evolution using the spherical collapse 
model. This simplification leads to the inaccurate con- 
clusion that a local bias at formation stays local. That 
gravitational evolution generates non-local bias can be 
seen from the results of [35 in second-order in perturba- 
tion theory, although this particular point was not noted 
in that work. The best known example is the limit of 
this result when the "formation time" is taken to be at 
the far past, the so-called local Lagrangian bias, and was 
first emphasized in [40] and further explored in |41j . In 
this paper we develop a formalism that contains all these 
results in particular limits. Moreover, it extends them 
i) to higher-order in perturbation theory, ii) to include 
non-conservation of tracers (arbitrary formation rate and 
merging), iii) to consider biased tracers that do not flow 
with the dark matter (velocity bias). Non-local bias is 
particularly interesting in view of the fact that the local 
biasing prescription does not seem to agree well with sim- 
ulations 42, 43]. Our model of the non-locality generated 
by evolution gives a well-motivated model for non-local 
bias. 

This paper is organized as follows. In Section [IT] we 
develop a formalism to generalize previous work on bias 
evolution to include velocity bias. We show that gravi- 
tational evolution induces a quadrupole, and hence non- 
locality of bias, on large scales. If velocity bias is present, 
then a dipole is also induced. We illustrate these effects 
for the case of the evolution of initially scale independent 
local bias. 



Section III shows that, when there is no velocity bias, 
then the same results can be obtained from a Lagrangian 
formalism, provided the initial conditions are treated self- 
consistently. In so doing, we demonstrate that Eulerian 
and Lagrangian treatments do, in fact, yield the same bis- 
pectrum; we discuss this in the context of what appear to 
be contradictory statements in the literature. Section [TV] 
studies bias evolution when comoving number densities 
are not conserved, either because of merging, or because 
of the formation of new objects. In this case also, no 
dipole contribution is generated if there is no velocity 
bias. In Section [V] we extend our calculation to third or- 
der (for the case of no velocity bias), and show that the 
structure of the non-local bias generated is most easily 
described by Galileon fields, with a dipole arising from 
breaking of Galilean invariance of the bias relation when 
there is velocity bias. Appendix [A] makes the connection 
between the conserved and non-conserved non-local bias 
in the most general terms. 



A comparison with simulations is done in Section VI 



where we use the results of previous sections to motivate 
a search for correlations between the halo overdensity 
at fixed matter overdensity with the different non-local 
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fields that our calculations singled out, finding signatures 
of non-local bias and its dependence on halo mass. In 
Section |VII| we discuss the impact of non-local bias on 
the bispectrum, and quantify the magnitude of non-local 
bias in simulations. A final section summarizes our con- 
clusions. 

Where necessary, we assume a flat ACDM cosmology 
with f2 m = 0.25 and J7a = 0.75. In this paper we 
use galaxies, halos, and biased tracers interchangeably. 
Those readers interested in skipping the technical details 
and focussing on the main results, the detection of non- 
local bias in the simulations and their implications, can 
jump directly to Section |VI[ where the main results de- 
rived previously are summarized. 



II. NON-LOCAL BIAS GENERATION WITH 
CONSERVED TRACERS 

A. Conserved Tracers with Velocity Bias 

We start by generalizing previous results [351 133 on 
the evolution of a tracer density perturbation (galaxies 
or halos), under the assumption that they form at a sin- 
gle instant in time with local bias, and thereafter evolve 
conserving their comoving number density (we relax this 
assumption in Section IV I. In particular, we include the 



possibility that these tracers do not flow with the dark 
matter, and therefore have their own velocity field. To 
fully specify the evolution of their velocity field however 
one needs to make some assumptions, here we will as- 
sume that the tracers are massless so we can ignore their 
contribution to the gravitational potential which is only 
sourced by the dark matter. This is a reasonable ap- 
proximation for galaxies, since only about 20% of the 



matter density is in baryons and an even smaller frac- 
tion of baryons is in galaxies [44] . At the large-scales of 
most interest, we can neglect dynamical friction and any 
pressure contribution, so we effectively treat the tracers 
as a pressure-less ideal fluid moving under the gravita- 
tional force generated by matter perturbations. In many 
respects, our approach is very similar to the perturba- 
tion theory treatment of two-fluids in |45j . a connection 
we will make more explicitly below (see also [SHH]). I n 
section fVl] we will apply our results to dark matter halos 
in simulations. In this case we are effectively assuming 
that halos may be treated as test particles (represented 
by their center of mass) whic move in the gravitational 
field due to the full matter distribution (i.e. all other 
halos) . 

In what follows, we will make heavy use of results 
from perturbation theory (PT, see |48j for a review). 
See section [V] for a simpler approach (in real instead of 
Fourier space) that neglects velocity bias, but which goes 
to third-order in PT instead of the second-order calcula- 
tions we do here. We assume that our tracers (which we 
will henceforth denote as galaxies) are formed at a single 
instant, with a spatial distribution that is a local function 
of overdensity S, and a velocity bias that is linear. We 
thus have two density and two velocity fields, one each for 
matter and tracers, and equations of motion that follow 
from imposing conservation of mass and tracers (we go 
beyond conserved tracers in section IV I and momentum 



conservation describing motion under the gravitational 
potential that is sourced by matter perturbations. 

For a single-component fluid, mass and momentum 
conservation can be combined into a single equation for 
a two-component "vector" which simplifies obtaining the 
evolution of density and velocity fields at once [I§1 [50] . 
In what follows, we generalize this to a four-component 
vector equation for our two-component model. That is, 
we consider 
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dr 



or 2 



d 3 fci d 3 fc 2 S u (k - k 12 ) a(k!,k 2 ) 0(kO 6(k 2 ), 
Hfi + -H'il m <) = - I d 3 fcid 3 fc 2 <5 D (k-k 12 )/3(k 1 ,k 2 )0(k 1 )0(k 2 ), 
d 3 fc! d 3 fc 2 J D (k - kia) o(ki,k 2 ) g (ki) <S g (k 2 ), 
d 3 fci d 3 fc 2 «y D (k - k 12 ) /?(k x , k 2 ) g (ki) g (k 2 ), 



(1) 
(2) 
(3) 
(4) 



where S and 9 are the density contrast and velocity di- coupling kernels a and j3 are defined as 
vergence of dark matter and S g and 9 g are the corre- 
sponding quantities for galaxies, r is conformal time, 
% = dlna/dr, and k i2 denotes ki + k 2 . The modc- 



a(ki,k 2 ) 



ki2 ■ ki 



/3(ki,k 2 ) 



fc? 2 (ki-k 2 ) 
2.- k-^ 



(5) 



We then introduce y — In D as the time variable, where 
D is the linear growth factor for the matter perturbations 
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satisfying 



d 2 D 
"dr 2 " 



H- 



dD 



-H 2 n m D = 0. 



(6) 



Since f 2 = f2 m , with / = dy/ din a, is a very good ap- 
proximation throughout the evolution |51j , the equations 
of motion Eqs. ([l]|4]) can be written in compact form by 
defining a four-component "vector" \I> as 



* = 




(7) 



linear velocity bias b* . Our choice of = means that 
we can restore the more general time dependence by re- 
placing y — > \nD/D* in all formulas below. The initial 
conditions can be handled conveniently by Laplace trans- 
forms. Taking the Laplace transform with respect to y, 
Eq. ([8| becomes 

CJ*a(w) - <j> a + tt ab y b (uj) = labc [* 6 * * C ]M, (12) 

where ^(a;) represents the Laplace transform of ^(y), 
4>a = ^a(y* = 0) is the initial condition and 



1 



x-\-ioc 



dw'f 6 (w')*c(w - w'), (13) 



which yields 

d y v a (k) + n ab v b (k) = labc ni^M, (8) 

where integration over and k 2 is implied and the en- 
tries of j abc are zero except for 



7121 — 7343 
7222 = 7444 



<5 D (k- ki 2 )a(ki,k 2 ), 
<5 D (k-k 12 )/3(k 1 ,k 2 ), 



(9) 
(10) 



and the matrix fl ab reads 
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(11) 



In this section we assume that galaxies are formed at 
a single epoch y* = with linear density bias b\ and 



for some x in the region of convergence of Collecting 
the terms linear in ^, we have 

^ a (uj) = a ab (cj)((f> b + 7f, cd * c * *d(w)), (14) 

with a ab — (ujI+Q)^ . Finally, we go back to the y-space 
by the taking the inverse Laplace transform 



^ a (y) = g a b(y)(j>b 



where g a b(y), called the linear propagator 
given by 



dy'g ab {y - y')ibcd^c{y')^d{y'), 

(15) 

ED, is 



Ja b (y) = 



l 

2iri 



duj a ab (uj)e u 



(16) 



where £ is a real number larger than the real parts of the 
poles of a. We then have, 



/ | e - 3a /2 + | e y 



9ab 



-1 



5 , 5 

. 3g-3y/2 _|_ 3_ e y 

I fe" 3 ^ 2 + 

3 e -3y/2 , 3e^ 
5 e ^ 5 



-3y/2 



| 5 e -3y/2 + 2gj/ 







(-2 - fe" 3 ^ 2 + 2c-y/ 2 + *f) 1 2(1 -c-y/ 2 ) 

3 e -3y/2 _ c -y/2 + q Q -y/2 



(17) 



We note that the 2x2 block in the upper left corner is 
the same as the linear propagator for dark matter derived 
in [49j [50] . The linear propagator satisfies the relation 

g a b(yi + V2) = g ac (yi) g C b(v2), (18) 

which is the expected law for linear evolution generalized 
for arbitrary mixing of growing and decaying modes. The 
linear propagator has the usual (matter only) growing 
and decaying modes, 




i.e. g a b{y)e b 1) = e y e£' and g ab (y) e [ h 2) = e 3y / 2 e a 2) . In 

(3) 

addition, there is an iso-density decaying mode e a and 
an iso-density- velocity decaying mode e a 4 \ 
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-W 2 \ 




1 


2^2 \ 


P (3) _ 


\ 












0/ 




V 


LJl } 



(20) 



where here we restored (following |45| ) temporarily a con- 
tribution to the overall mass density fraction of lj\ for 
matter and w 2 for galaxies + cj 2 = 1). The first 

eigenmode here satisfies g a b(y) ej, 3 ' = el , correspond- 
ing to a constant mode with zero total density perturba- 
tion, while the second eigenmode satisfies g a b(y) e^ 4 "* = 
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gT v I 2 and corresponds to a vanishing total density 
and total velocity divergence perturbation. Our assump- 
tion of tracers as test particles (massless) means we have 
set uj\ — 1 and cj 2 = in our approximation. In the more 
general case, the same results we find here will apply with 
small corrections proportional to w 2 (see [45]). Note that 



(1 2) 

the standard eigenmodes e a ' are of course independent 
of LOi as they correspond to in-phase motion of the two 
fluids as if they were one. 

We have transformed the equations of motion Eq. ^ 
into an integral equation Eq. ( |15| with explicit depen- 
dence on initial conditions <p a that can be solved pertur- 
batively. To specify the initial conditions we assume that 
they can be expanded as follows, 



a (k)=]T / d 3 qi ... d 3 g „5 D (k-q 12 ...„)Z(")( qi ,...,q n ) <J (qi)...<So(qn), (21) 
n J 

where 5q is the initial dark matter density contrast, and the vector \E' a (k,y) can be similarly expanded as 

* (k,y)=£; J d 3 gi ...d 3 ( 7n( 5 D (k-q 12 ... n )/Ci")(q 1 ,...,q n , 2 /) ( 5o(qi)-.-'5o(qr 1 ). (22) 



the K.^ kernels, 



Putting Eqs. (21 1 and (22 1 into Eq. ( 15 I, and collecting terms of the same order in So, we get a recursion relation for 



rv 

£i n) (qi> ■■•,<bi,y) = 9ab(y) T-i n) (qi, . . . , q„) + V / dy'g ab (y - y') 7 6cd (ki = qi...j,k 2 = q J+ i. 

3=1 Jo 

x IC ( C J) (qi , . . . , q^ , y') JC { ^~ 3) (q J+ i , . . . , q n , y'). 



(23) 



Note that the kernels K.^ obtained from Eq. (|23|) arc 



not symmetric in the arguments q^, but they can be 
symmetrized afterwards. Only the symmetric part con- 
tributes to 



B. Generation of Non-Local Bias 



"Initial Conditions" at Formation 



We now explore the results of Eq. (23) to study the 



generation of non-local bias by gravitational evolution 
from local-bias initial conditions. That is, we assume 
that biased tracers at formation time t* (corresponding 
to growth factor D* and y = 0) can be written as a local 
function of matter density that is then Taylor expanded, 



_k S k 



b* 0* 



(24) 



where we assumed that the tracers have velocities that 
are only linearly biased with respect to matter. In the 
examples below we assume 6 V = const, but our results in 
this section also apply if the velocity bias is /c-dependent. 
Small-scale velocity bias has been seen in simulations 
[5"2Tf5"4"] . at the 10% level. At large scales is predicted 
by peak theory [55l [56] although in a statistical sense, 
i.e. peaks move locally with the dark matter but their 
statistics can be thought of as if there is a velocity bias 
that is fc-dependent and goes to unity at very large scales 



as k 2 . But the situation for baryons, and therefore galax- 
ies (as opposed to halos), can be somewhat different, e.g. 
at early times z > 100 the relative velocity between the 
dark matter and baryons is typically supersonic |57l 158) , 
and there might be a non-trivial component to the re- 
lation between dark matter and baryons on large scales 
due to isocurvature modes, see e.g. [59] . 

We assume that the matter is in the growing mode, so, 
to linear order, the initial conditions kernel at formation 
time is given by, 



4 X) - I 6 j J = 4 X) + [(65 - 1) + 2(6; - 1)] 4 3) 

+ (6;-l)4 4) . (25) 

Here we have expanded the initial conditions in terms 
of the eigenmodes of the linear propagator (recall that 
u)x = 1 and ll>2 = in Eq. 



20). This makes clear that 



density bias excites the iso-density decaying mode and 
that velocity bias excites, in addition, the iso-density- 
velocity decaying mode. At second-order we have, 



T( 2 ) - 



/ F 2 (ki,k 2 ) 
G 2 (ki,k 2 ) 
b* 1 F 2 {k 1 ,k 2 )+b* 2 /2 

V &*G 2 (k!,k 2 ) 



(26) 



where F 2 and G 2 are the second-order kernels that de- 
scribe the density and velocity divergence to quadratic 
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order in the linear matter fluctuations (see Eqs. 34p6 
below for the multipole expansion of i 7 ^)- They are gen- 



crated by the second term in Eq. ( 23 ) during time evolu- 
tion up to y = for matter fluctuations (which satisfies 
the same equations for y < as Eq. ( 23 1 restricted to 



a = 1,2). They are functions of the wavectors through 
Ibcd because gravitational evolution is nonlocal. 

At time t > t* (or redshift z < z*) S g will be, pre- 



cisely because of ^bcd in the second term in Eq. ( 23 ) , a 
non-local function of 5. We therefore are interested in 
separating out the non-local contribution to the galaxy 
density perturbations, 



rNloc _ 



Hocal 



(27) 




where <5g° cal is a local function of 5 and thus can be ex- 



pressed in terms of local bias parameters (the evolved 
version of the b* k 's) . We will in fact construct 5g loc and 
<5g° cal order by order, e.g. by first substracting linear 
bias and analyzing local and non-local contributions at 
second order, then substracting local quadratic bias and 
analyzing what happens at third order, and so on. 

In this section we take the first step in this analysis: 
to quantify the non-local contributions to second order. 
To do so, it is convenient to define the field 



X = S g - b 1 S 



(28) 



which neglects all k > 1 terms in S g loc of Eq. ( 27 1 . We 
can then study x a ^ second order and decompose it in 
Legendre polynomials Ve(p), 



X 



(2) 



(x) = J d 3 fc 1 d 3 fc 2 e- lk - x ( 5 (k 1 )<5o(k 2 ) 

OO 



(29) 



1=0 



where fi = k\ ■ k%. A local contribution to x^ should 
be proportional to <5q, and therefore corresponds to a 
monopole (£ = 0) contribution with xo independent of 
ki. Any £ > piece cannot be written as local functions 
of 6 and thus will be entirely the result of non-local con- 
tributions. As we shall see, a quadrupole contribution 
(£ = 2) is inevitable for biased tracers, and velocity bias 
generates in addition a dipole (£ — 1). We shall not go 
beyond second order here, see Section [V] for results to 
third order (see also Appendix |A|) , and the next step in 
the construction of Eq. (27). 



2. Evolution of Linear Density and Velocity Bias 



We now turn to the solution of Eq. (23). To linear 
order, only its first term contributes. From the decom- 



position into eigenmodes of the propagator, Eq. (25), we 
can read off the evolution of each field, which is precisely 
of the local form given in the initial conditions but with 



FIG. 1. The evolution of linear density bias 6i (top three 
lines) and velocity bias b v (bottom lines) as a function of the 
scale factor a with the initial values bl — 2 and b * = 1.1, 1 
and 0.9 respectively. A velocity bias larger than 1 slows down 
the decay of density bias slightly, while velocity bias less than 
1 speeds it up. 



a prescribed time dependence. For matter density and 
velocity fields we have linear growing-mode evolution, 



4 1) (y)=4 1) (y) 



(30) 



while, for density and velocity bias we have, respectively, 

h ee = i + (&*- i) e -* + 2(6; - i) c -»(i -c-y/ 2 ), 



^\y) 



(31) 



and 



b v ^^fl = l + (b* v -l) c-*y/\ (32) 

K-2 \V) 



Note that when there is no velocity bias, b* 



recover from Eq. (31) the well-known result 

(bt 



l)e" 



1, we 



(33) 



that density bias asymptotes to unity in the long-time 
limit if the comoving number density of tracers is con- 
served. On the other hand, our generalization to b v =/= 1 
does not agree with recent assumptions about the evo- 
lution of peaks in the initial density field [56] . which do 
not show the presence of the iso-density-velocity m ode 
contribution that gives the e _3y / 2 decay in Eqs. (31 
32 ) . This disagreement results from different assump- 



tions. Peaks move locally with the dark matter but their 
velocity statistics can be thought of as if they had a statis- 
tical velocity bias that remains constant with evolution. 
Because of this difference in treatment, the peaks calcu- 
lation cannot be directly compared to what we do here, 
although it is important to clarify which treatment is a 
more accurate description of velocity statistics of tracers. 
We hope to report on this in the near future. 
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Figure [T] shows b\ and b v as a function of the scale 
factor a. We have set b* = 2 and 6* = 1.1,1, and 0.9. 
Note that 6* > 1 slows down the relaxation of the density 
bias slightly while 6* < 1 speeds it up. Velocity bias also 
relaxes to unity eventually. 



3. Quadratic Order: Emergence of Non-Local Bias 

Because the vertex ^bcd is quadratic in ki • k 2 , only 

I < 2 multipole moments will be present. For the matter 

(2) 

density field, the multipole expansion of K\ (y) reads, 
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(2) _ - p 2y 

21 

1 (ki 

2 U 2 



rw _ _ „ 
^1,^=0 ~~ ~ 13 



K (2) 
^11=1 



K, {2) - — c 

/v l/=2 — 21 



2y 



(34) 
(35) 
(36) 



These correspond to the multipole expansion of 
e 2y F2 (ki , k 2 ) . The monopole represents the second-order 
nonlinear growth in the spherical collapse dynamics, the 
dipole the transport of matter by the velocity field, and 
the quadrupole describes tidal gravitational effects. 

The multipole moments for the galaxy density pertur- 
bation to second-order, \\ ■> are given by 



(2) 

Xo 



(2) 

Xi 



21 65 + 2T £v 



jfel / 
k 2 h , 
4 4 



3 + Uc y/2 - 14e, 
1 + 2e y/2 



^ 2) = -2T £5+ 2T^ 



ef - 1 
12 + Ue y/2 + 7e 



21e^ - 
ev - 1- 
(37) 

(38) 
(39) 



where e$ and e v are proportional to density and velocity 
bias, respectively: 



e 5 = (h-l) e y (e y -l), 



(6v-l) e» (e^ 2 -l), 



(40) 



(41) 



and vanish for fully unbiased tracers. They are also de- 
fined to be zero at formation time (y = 0), leaving only 
the prescribed monopole from local bias. In the long-time 
limit (y — > 00) they asymptote to 

e s ^e y [(b* 1 -l)+2(b* v -l)}, e v ^(b*-l). (42) 

Thus, the effects proportional to dominate, but they 
are suppressed by another factor of c y compared to the 
usual second-order effects. Finally, note that when there 
is no velocity bias (e v = 0) Eqs. ( 37p9 ) reduce to, 



(2) 

Xo 



21 



£5, 



(2) 
Xi 



0, X2' 



(2) 



21 



e<5, (43) 
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quadrupole, b ¥ ' =1 
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FIG. 2. Emergence of non-local bias from local-bias initial 
conditions, as quantified by the ev olution of the ratio of the 
galaxy multipoles Xe 

r> (Eqs. |37|39| to the corresponding mat- 



34p6 l 



At formation (a = 0.2), bias is 
0.5, i.e. there is only a monopole at 



ter multipoles (Eqs 
local with bl = 2 and~6j 
second-order. However, a quadrupole (three bottom lines) is 
generated at later times. If there is velocity bias, then a dipole 
is also generated (three middle lines) . The three lines for each 
multipole correspond to different choices for the initial veloc- 
ity bias: b* = 1.1 (solid), 1 (dashed) and 0.9 (dotted). 



with es = {bl — l)(e y — 1). In this case there is no dipole, 
and the induced quadrupole and monopole are opposite 
in sign. The induced structure when there is no velocity 
bias for conserved tracers is further explored in Section 
|V| to third-order in PT, and in Appendix [5] for the non- 
conserved case. 

Figure 2 shows the evolution of the multipoles x±}_ 
(normalized by their dark matter counterparts, Eqs. |34[ 
[36| as a function of y for three different choices of velocity 
bias, b* = 1.1 (solid), 1 (dashed) and 0.9 (dotted). We 
see that even though the bias at formation (a = 0.2) is 
local (only a monopole is present), higher-order multi- 
poles get generated. If there is no velocity bias then only 
a quadrupole gets generated; if b v ^ 1 then a dipole is 
also generated (with sign determined by b* — 1). All of 
these normalized multipoles eventually relax to zero be- 
cause the galaxy multipoles grow more slowly than those 
of the dark matter. If b* = 1.1, the relaxation of the 
monopole is slowed down, whereas the quadrupole re- 
laxes faster; the opposite holds for 6* = 0.9. We see 
that even for a significant velocity bias of 10%, the gen- 
erated dipole is only 10% of that in the dark matter. A 
dipole contribution in galaxy bias can enhance the shift 
of the BAO peak in the correlation function [SO], but 
since the dark matter dipole effect is at the percent level, 
velocity bias is unlikely to change this in any significant 
way except possibly for the very highly biased tracers. 
See Section I VI I for more discussion on the effects of such 
dipole term from numerical simulations. 

Thus, we see that non-local bias is inevitably in- 



duced by gravitational evolution, and that the locality 
assumption cannot be self-consistent. In practice, be- 
cause galaxy formation happens in a continuous fashion, 
we don't expect locality to be valid at any time, even if 
the formation bias were local. We explore this in Sec- 
tion [iVj In addition, there is no reason to expect the bias 
at formation to be purely local, even for dark matter ha- 
los, since the barrier for collapse is known to depend on 
quantities other than the overdensity (e.g. [35]). 



III. COMPARISON WITH LOCAL 
LAGRANGIAN BIAS 

In this section, we would like to compare our results 
with those known from the literature on local Lagrangian 
bias, which can be thought of as a particular limit of our 
results when formation time is at the far past (z — > oo) 
and there is no velocity bias. While such calculations are 
usually done in Lagrangian PT (see e.g. 001111]), clearly 
one should obtain the same results if done in Eulerian PT 
as we have used so far. It is however instructive to redo 
this calculation in a Lagrangian description and compare. 

Since there is no velocity bias and tracers are con- 
served, the continuity equations Eq. ([I]) and Eq. ^ can 
be used to relate <5 g to S through the matter velocity di- 
vergence field HO], 

dHi + 5) = y u= dHi + s g ) (44) 



(It 



(It 



where we used the Lagrangian or total derivative follow- 
ing the motion of a fluid element, 



d _ d 
dr dr 



u - V. 



Upon integration of Eq. ( 44 1 , we get 



ln[l + <J g (x)] = ln[l + <J(x)] + /(q), 



(45) 



(46) 



where /(q) is a function depending on initial fields at 
the Lagrangian coordinate q related to the Eulerian x 
through the displacement field 



x(q,t) =q + *(q,t). 



(47) 



In terms of the initial condition /(q) is clearly given by 
ln{ [1 + 5 g (q)] / [1 + S(q)] } and thus we have, 

where the Lagrangian fields are evaluated at the ini- 
tial time £*. This is the same result as that given 
in [ini Ell HEj except for the denominator (1 + <5(q)), 



which was implicitly dropped in those works (it is how- 
ever included in |61|). However, to reproduce the decay- 
ing modes found in the previous section, this denomina- 
tor is required. 

For comparison with the results in the previous sec- 
tion, we now assume local Lagrangian bias in Eq. ( 48 ) to 
quadratic order, 



l + J g (x) = [l + *(x)] 



l + bfS(q) + (bf/2)5 2 (q 
1 + S(q) 



(49) 



where the Lagrangian bias parameters bf are the equiv- 
alent to the parameters b* in the previous section. To 
linear order, we can assume q ~ x in this equation, but 
to go to second order we need to include the displace- 
ment field to first order, i.e. in the Zel'dovich approxi- 
mation (hereafter ZA, [55]). This is given by 



*(q,t) 



D(t) - D{U) 
D(U) 



■ ' p ) ^ 



a ik q 



(50) 

where we have used the fact that, in the ZA, the de- 
caying mode is constant [49]. Note that ^(q, t*) = at 
formation time t*, as it should. We emphasize again that 
this decaying mode, which is often neglected in the liter- 
ature, must be included if one wishes to fully reproduce 
the results in the previous section to second order. 

To linear order Eq. Q reads ~ S^+bf S^-SP, 
where all fields have the same argument x. Therefore we 
deduce the Eulerian linear bias bf 



bf = l 



bf - 1 



(51) 



which is Eq. (33 1, with y = hi(D/D*). This seems differ- 



ent a priori from the often quoted relationship between 
linear Eulerian and Lagrangian bias bf = 1 + bf. The 
reason is twofold: first, it is customary to define the La- 
grangian bias with respect to the extrapolated linear den- 
sity field 5^ rather than the Lagrangian density field 5* 
as we have done here, so the more standard definition is 
instead 



bf^bf 



~D 



(52) 



and second, if we neglect the third term in Eq. ( 51 ) com- 
ing from the denominator in Eq. ( 48 ) , then we recover 
the familiar bf = 1 + 6^. This second step is justified 
for objects that are not arbitrarily close to unbiased in 
which case as t* — » and D/D* — > oo, bf increases with- 
out bound for objects with fixed bf, so that bf ^> 1 in 
this limit. Although this step is unjustified for unbiased 
objects for which bf = 1, keeping only this term does no 
harm as its contribution to Eq. (51 1 vanishes as t* — > 0. 
To second order in PT, Eq. (49) gives, after using 



Eq. (47) 
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4 2) w = 



i 



b\ - 1 



8<V (x) 



b\ - 1 
D/D* 



V<5 (1) (x) 



(«5«(x)) 2 



which upon Fourier transform, after using Eq. (50), we can write as the quadratic kernel for galaxies 

2 



~D 



K. 



(2) 



5 b^-1 
7 + (£>/£>*) 



1) 



1 



6^ - 1 



2 U 2 



k2 



i 



6^-1 



(53) 



(54) 



wher e we h ave used the second-order matter results 
Eqs. (34 36 1 , and /i = ki • k 2 . It is easy to check that this 
equation agrees with Eq. (43) with e$ = (bf — l)(e y 



1) 



after using that x^ — ^-3 — bf K\ with bf given by 
Eq. J5l]) and y = ln(D/D*). 



We can now take the limit t* — > in Eq. (54) as for 



the linear result above and compare with the results in 
the literature. Now we need to redefine the quadratic La- 
grangian bias in terms of the present density fluctuations, 
as done for the linear bias in Eq. (55), 



b\ = b 



~D 



(55) 



and assuming bf ^> 1 as before we then get for Eq. ( 54 ) , 




2 \h 



7 



M 2 , (56) 



which agrees with Eq. (8) in [41]. Note however that 
in [41] it is argued that the dipole term (proportional 
to b\ in Eq. 56) is a new feature of local Lagrangian 
bias as opposed to local Eulerian bias. This interpre- 
tation is not correct: in local Eulerian bias the second- 
order galaxy kernel is, apart from local contributions of 
quadratic bias L that of the matter multiplied by linear 
bias bf = 1 + bf, so the precise amplitude of the dipole 
agrees with that in Eq. ( 56 1 . Subtracting this local Eu- 
lerian piece to construct x^ = ~~ bf K-i^ 1 shows 
that the new contributions are indeed of the form given 
by Eq. (43). That is, the new qualitative contribution 
is a quadrupole term, not a dipole. As we showed in 
the previous section, an additional dipole will only ap- 
pear if there is velocity bias. The physical reason for this 
(breaking of the Galilean invariance of the bias relation) 
is discussed in Section [V] 

Thus, we have shown that, in the appropriate limit, we 
reproduce known local Lagrangian bias results. However, 
our approach in the previous section is more flexible as 
it does not require "formation" to be in the distant past, 
and it also allows for velocity bias. On the other hand, 
we have, so far, assumed that the comoving density of 
tracers is conserved. We now discuss how to go beyond 
this assumption. 



IV. NON-LOCAL BIAS GENERATION WITH 
NON-CONSERVED TRACERS 

Galaxies form at a range of redshifts and merge. So it is 
important to extend the previous results to the more real- 
istic case when the comoving number density of galaxies 
changes with redshift due to some arbitrary source field 
j, which effectively includes the effects of galaxy forma- 
tion and merging. Our description here is similar to |37j 
(see also [S3] IM]), but we shall extend the analysis to 
higher order in PT. For simplicity here we assume that 
the bias at formation is local (as we have done so far), 
Appendix [A - ] discusses what happens in the most general 



case (see also Eq. 114 below). The evolution equation for 



the physical galaxy number density n g now becomes 



& 

dr 



3Hn g + V • (n g u) = AUj{p). 



(57) 



Note that we factorized the source term into two com- 
ponents, A and j, where A roughly parametrizes the 
epoch of galaxy formation (e.g. following star forma- 
tion history) and j describes the effects of dark matter 
on the formation and merging of galaxies. Nonetheless, 
we stress that our main results are independent of the 
detailed functional form of A and j, and we use the as- 
sumed functional forms only to make the plots shown in 
Fig. [3] For example, A(t) can be a log-normal profile 



1 _(l»«-l»«o) 

A{t) = — e 
a 6 



(58) 



where a and 00 are free parameters. For j(p) we take a 
simple quadratic form 



j(p) = on — + "2 f — 
Po \Po 



(59) 



where a\ and a-i are model parameters, and pa is the 
average matter density today. Appendix \K\ considers the 
implications of j depending on non-local functions of 8, 
or other fields. 



The second term in Eq. ( 57 ) can be eliminated if we use 



comoving rather than physical number densities, rig 
a 3 n„, so we have 



(c) _ 



dn, 



(c) 



1 



d In a H 



+ -\7-(n^u)=A^j(p), 



(60) 
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FIG. 3. Local bias parameters at formation bl (top-left panel) and b% (top-center), comoving number density (normalized by 
present value, top-right), linear bias (bottom-left), second-order galaxy bias quadrupole X2 an d monopole Xo 2 ' (normalized 
by dark matter values, bottom center and right panels) as a function of scale factor a. Each panel shows three sets of values of 
ai and CK2, corresponding to {01,0:2} = {4, 1} (solid), {1, 1} (dashed), and {1,4} (dot-dashed). To describe galaxy formation, 
we have used the toy model in Eq. (58 1 with 00 — 0.2 and characteristic galaxy formation time do equal to 0.3 (blue), 0.5 (red) 
and 0.7 (green) respectively. 



where A^ c > = a 3 A. We then write, 



(61) 



and solve Eq. ( 60 ) by perturbation theory. We will also 



assume that there is no velocity bias, so that galaxies 
and matter share the same velocity field, which is known 
from solving the evolution of matter. We then expand 
the source term on the right hand side of Eq. (60) to 
second order 



1 



(62) 



This invites us to interpret j'(p)p/j{p) and j" (p)p 2 /j(p) 
as the instantaneous formation bias b\ (t) and b* 2 (t) of the 
galaxies formed (or destroyed) at time In a 



mi 



W) 



j"(p)p 2 
m 



(63) 



For example, the form of j in Eq. (59 1 gives 



b\ (a) = 1 + 



«2 



«ia 3 + a-i ' 



2a 2 



a 2 



(64) 



which imply a simple relation b\ — 1 = 6^/2. This relation 
only holds at the formation time, as we will see evolution 
inevitably generates non- locality and breaks this. Note 
that in Eq. ( 64 ) the bias parameters have a pole when 
0L2 is negative. This means that quadratic approxima- 
tion is no longer valid and the higher order terms in the 
expansion are important. 

Figure [3] shows in the top left and center panels b\ and 
b* 2 for three sets of values of ai and corresponding 
to {ai,a 2 } = {4,1} (solid), {1,1} (dashed), and {1,4} 
(dot-dashed). As we mentioned before, these choices are 
just illustrative with no special physical significance, but 
serves to show a range of possibilities for the local biases 
at formation. 
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A. Background Solution 



We now look for the evolution of the galaxy comoving 

- (c) 

number density rig by solving the background equation, 



which gives for the effective linear bias, b\ = D g /D 
1 /"™e c) 

h(t) = l + -M- dn, (bl-l)D„ (72) 

nh D Jo 



dn g c) 
din a 



(65) 



where p is the mean matter density. The solution is 



In a 

n\ = I d(\na)A^j(p), (66) 

'lnOini 



iW = 



where we have assumed that there are no galaxies at 
t = tjxa. The top right panel in Fig. [3] shows the resulting 

-(c) 

n g for three sets of parameters a% and a.2 in three differ- 
ent backgound profiles A. As expected the background 

_(c). 

number density n g is predominantly determined by the 
profile A. 



B. First-Order: The Evolution of Linear Bias 



where quantities inside the integral are evaluated at 
a time when the comoving number density of galaxies 
equals . This equation is the generalization of Eq. ( 33 ) 
(with y = In D / 'D„ ) for when galaxies form during a 
broad range of redshifts and are not necessarily con- 
served. We have assumed for simplicity that the rela- 
tionship fig '* to In a is one-to-one, if not (because merg- 
ing overcomes formation at some periods of time), one 
should just sum up over all contributions at a given value 
of n in Eq. ( 72 1 . The key idea is that the simple result 



of the conserved-tracers case, Eq. (33 1, gets generalized 



by simply weighting by instantaneous comoving number 
density. 



The form of Eq. ( 72 ) is not practical for observational 



To first order in PT we write n 



(c) 



4 C) (1 + 4 1; ) in 



Eq. ( 60 ) , and using the background evolution in Eq. ( 65 1 , rewrite Eq. ( 72 ) as, 



purposes since it is not feasible to trace back all the way 
to when there were no galaxies. A more useful form is to 
use initial data at some (high) redshift when the growth 

factor was Di and there was some non-zero comoving 

_( c ) 

number density n ■ with linear bias b\i . One can then 



we get 
dsP A^ 



;(c) 



d In a 



bi{t) = 1 + 



Di 



nl c) D 



(67) 

which using linear theory evolution for matter with 
growth factor D can be rewritten as 



dn* (&J-1) £>*, 
which follows from combining Eq. (72 1 at times ij and t. 



OS. 



(i) 



A( 



This also makes rather clear that if the comoving number 
density does not change between ti and t (n 
the last term vanishes and one recovers Eq. 



(c) 



-| —j(p)$W = (/ -| ^xf(p)p)D Sq. (68) tracers are not conserved, we see that the second term in 



9 In a ' fgyv « ^ fl (c) 



= 4 C) ), 

. When 



Eq. ( 73 ) is simply the usual (conserved case) decay of bias 



where / = d In Dj din a. Looking for a solution of the 
form 8g = D g (t) Sq, where D g is the linear growth factor 



modulated by the evolution of n e with the result that 

this decay is slowed down by merging (decrease of n g °^). 
In addition, the third term provides an extra contribution 
that depends on the bias of galaxies that are formed or 
lost to merging, and the sign of this contribution depends 
on whether formation or merging dominates. 

Another way to see these effects is to go back to 
Eq. ( 69 ) and write an explicit differential equation for 



for galaxies, and after using Eq. (65), we arrive at 
d 



(n^D g ) = (n^f + A^f(p)p)D. 



din a s 

And since n g — initially, we have 



(69) 



the linear bias, by replacing D g with b\(t) D, which gives 



1 



D g (t) = D(t) + ~ I d(\na)A^D(-j(p)+f(p)p), 

fl g Jlna ini 

..(70) 

which agrees with [37 after some simple redefinitions. 
We can express Eq. ( 70 ) in a more physical way that 
makes clear the independence of the detailed form of 
and j, by using the comoving number density of galaxies 
as the integration variable, 



d 

d In a 



(h - 1) = -(h - l)f 



din n L 
din a 



(c) 



(h-bt). (74) 



The first term on the right hand side is the usual term 
that drives the decay of bias due to the growth of the 
large scale structure. Indeed, neglecting the second 
source term one can solve Eq. ( 74 ) to recover Eq. ( 33 1 . 



D g (t) - D(t) + 



1 

-(c) , 



dn, (b\ - 1) D, 



(71) 



The second source term may speed up or slow down the 
decay of bias compared to the conserved model, depend- 
ing on the signs of (dn^/da) (i.e. whether formation or 
merging dominates) and (bi — b* ) specifying whether the 
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galaxies created or destroyed are more or less biased than 
the overall bias. 

In the bottom left panel of Fig. |3j we plot the the 
evolution of b\ for several choices of our parameters. It 
is interesting to note that as the characteristic galaxy 
formation time ao increases, the difference between the 
three models represented by {cti, a^} (different line types 
in Fig [3J becomes more marked. From Eq. (64) we see 



that if a2 is small compared to ai, then there is some 
epoch where (61 — b*) is positive, and this can accelerate 
the decay of bias. For the model {ai, a 2 } = {4, 1} (solid 
line in Fig. [3| , (61 — b*) becomes positive roughly when 
a > 0.5 and the decay of bias is speeded up. This is 
particularly apparent for the model with ao = 0.7 (solid 
green) as galaxy formation occurs when (bi — b\) is pos- 
itive. 

Current galaxy surveys are deep enough to test the evo- 



lution of linear bias, something that will be done much 
more precisely in the near future. Comparison with the- 
oretical predictions are often limited to the conserved- 
tracers case, Eq. (33), typically showing a decay of bias 
that is faster than predicted by this formula [6"STl71| . 
These deviations are interpreted in terms of merging 
and/or disruption of galaxies, and our more general for- 
mula, Eq. ( 73 ) , can help understand the implications of 



such measurements for merging/formation rates. 



Second-Order: Generation of Non-Local Bias 



To find the second-order solution, we write n 



(c) 



nl c) (l 



r(l) 



(67), we obtain 



in Eq. (60), and using Eq. (65) and 



(2) 



d In a 



+ ^(p)4 2) = ~v • u( 2 > - Iv • (4 1} u (1) ) + A> 



n 



(c) 



(75) 



(2) 

We then rewrite the source terms as before, in terms of solution of the equation x? = 0. That is, 

_(c) 

rig , b\ and b* 2 > substract the linear bia s contribution at 



second-order to define x as in Eq. ( 28 ) , and decompose 



it in terms of Legendre polynomials as in Eq. ( 29 ) . This 



(2) 

gives differential equations for the x\ ■ These can be 



D / dn,D,(bl - 1) 
Jo 



dn*D^-l). (79) 



integrated as before (see Eq. H by using n g c) as the Differentiating this equation with respect to n g c) leads to 



integrating factor, to obtain (compare with Eq. 43 1 



dn»D„(&J - 1) =0. 



(80) 



(2) 

x\ = 



21ni c) Jo 



dn*D*(D -D m )(bl- 1), (76) 



(2) 

Xi 

(2) 

Xo 



= 0, 



~x 2 



(2) , 1 



-(c) , 

ni Jo 



(77) 
(78) 



Note that the dipole moment is exactly zero, as in the 
conserved case, irrespective of the functional form of the 



source term. In fact, we can recover the results of Eq. ( 43 ) 
by inserting od(«* — n^^ng; in the integral and not- 
ing that = for the conserved-tracers case, and 



(2) 



1 in Section [ill Figure^ shows the evolution of Xi 



D * 

(bottom center panel) and \o (bottom right), normal- 
ized by the corresponding dark matter multipole coeffi- 
cients in Eq. ( |36| and p4| , respectively, as a function of 
scale factor a for the different choices of the parameters 
that describe A and j. 

It is interesting to ask if we can avoid inducing non- 
local bias by a proper choice of sources. Since locality is 
violated by the induced quadrupole term, we look for a 



Together with Eq. (71), this implies that D g = D, i.e. 
the galaxies must be unbiased. Hence biased tracers 
generically have non-local bias irrespective of their merg- 
ing/formation history. 

In Appendix[X]we study non-conserved tracers in more 
generality and show that the non-local bias can be related 
in a simple manner to that in the conserved case to ar- 
bitrary order. Since the lack of conservation of tracers 
does not qualitatively change the structure of the non- 
local bias, we now go back to the conserved case and 
discuss the basic features that determine this structure. 



V. STRUCTURE OF THE NON-LOCAL BIAS 
TO THIRD ORDER 

In this section, we derive the form of the non-local bias 
to third order. In principle, the general solution Eq. (23) 



developed in Section|II| (and generalized to non-conserved 
tracers in Appendix [A]j would allow us to write down the 
solution immediately. However, as we see in Eq. (43), 

(2) . . . . 

x\ exhibits an interesting structure when there is no 

velocity bias (preserved when going to non-conserved 
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tracers, see Eqs. 76 
space solution is rat. 



78). The full third-order Fourier- 
tier complicated, and naively apply- 
ing Eq. (23) to third order hardly gives us any insight. 



We find it more instructive to develop the solution in a 
different way. 

For simplicity we assume conserved tracers and no ve- 
locity bias, so the matter velocity field u controls the bias 
relation (see Appendix [A] for generalization beyond these 
assumptions). We assume the velocity field is potential, 
and, as discussed in Section [TTJ it is convenient to work 
with a normalized velocity field 



m 



v$ v , 



v 2 $ v 



(81) 



whose divergence agrees with density perturbations in 
linear theory, Oy — 8^\ We assume local bias at for- 
mation time t* (although see Eq. 114 below) 



h* h* 



(82) 



and look for the evolution of the bias relation, which can 
be divided into local bias evolution and induced non-local 
terms. Because the results to third order are compli- 
cated, we start by doing the calculation in the Zel'dovich 
approximation (ZA), where it is easier because the dy- 
namics of gravitational instability is local, and then ex- 
tend it to the exact dynamics. This serves to highlight 
the similarities and differences from the exact dynamics, 
and how its non-locality affects the bias relation. 

In the ZA, it is simple enough to see what to ex- 
pect. The dynamics is given by the displacement field 
if? that scales linearly with the growth factor (see Eqs. 47 
and 50 ) . The displacement field is related directly to the 
velocity potential by ^(q) = — V$ v , so the bias relation 
will be determined at large scales entirely by the second 
derivatives V,j$ v , which measure the variations in dis- 
placements that affect the clustering. The non-locality in 
the bias relation will then be determined by the scalar in- 
variants of Vij$ v , since the galaxy density perturbation 
is a scalar under 3D rotations and translations. 

In three dimensions, there are only three principal in- 
These are the "Galileons" [72] 



variants of 



e 2 (*v) 

£ 3 (<i>v) 



(V 2 $, 



v 2 $ v = ( 
(v, 3 $,) 2 

(V 2 $ v ) 3 

-3(V ij $ v ) 2 V 2 $ v 



2Vy$vVjfc$ v Vfci$v 



(83) 
(84) 

(85) 



Note that (£ 2 ($v)) = (^(*v)) = 0. On the other 
hand, one can construct a similar description through 
the invariants of the deformation tensor. In terms of the 
eigenvalues Aj of the deformation tensor 
these are 



h = Tr [V tj ] = Xi + A 2 + As, 
I2 = A1A2 + A2A3 + A3A1, 
J 3 = Det p tj \ = AiA 2 A 3 , 



(86) 
(87) 



making Q\ = I\, Q2 — ~2/2 and Q% = 6/3 to leading 
order (i.e. ignoring the distinction between q and x). In 

linear theory l\ = 8^ = but in general there is no 
simple relation between Ii and the Eulerian fields, unlike 
Q\. Since we work in Eulerian coordinates we will use 
the Galileons rather than the invariants of the deforma- 
tion tensor, although this distinction is not crucial. The 
invariants of the deformation tensor can also be related 
to other useful quantities in the theory of random Gaus- 
sian fields, such as the ellipticity and prolateness of the 
potential, see Appendix [B| 

Galileons arise in theories with the linear shift symme- 
try <E> V — > 4> v +V-x+C, which, in our 3D case corresponds 
precisely to Galilean invariance, since the equations of 
motion Eqs. (l|4) are Galilean invariant [73]. In theories 
of modified gravity, 4D Galileons arise as a remnant of a 
higher-dimensional Lorentz symmetry and Lagrangians 
containing Galileon interactions have equations of mo- 
tion that are precisely second-order [75]. In our case, the 
bias relation in the absence of velocity bias solely depends 
on $ v , which obeys second-order (in space) equations of 
motion. When there is velocity bias, however, there is 
a relative bulk flow between tracers and matter, which 
breaks the Galilean invariance of the bias relation, and 
this shows up as a dipole field at second-order in PT. This 
simple argument explains the basic structure of non-local 
bias. 

Therefore, we expect that in the ZA dynamics the bias 
relation when there is no velocity bias will have the form 
up to third-order, apart from local terms in 8 



rNlo 



= 72^2 ($ v )(l+/3<5)+ 7 3S 3 ($v) 



(89) 



with coefficients 7, oc (b\ — 1) and mixing j3 oc b\. Note 
that Q\ does not appear as 8 is already part of the local 
bias description, and the difference between Q\ = 8 V and 
8 can be written in terms of Galileons (see e.g. Eq. [94] 
below). When we extend these results to the exact dy- 
namics we shall see that new contributions appear that 
are nonlocal in $ v . The second-order term in Eq. (89) 



was already shown in Section [II] so the main task here is 
to show what happens to third order. 

To proceed, we start by substracting the local linear 
term, and writing an equation of motion for x = 8 g — b%8, 
which is simply (after using Eq. |33[ ) 

, (90) 



X = (6i-l)A + V-(xv), A = S- 



where a dot denotes a derivative with respect toy = In D 
The equation of motion for A is, after using Eqs. ( l][2) 



A+-eA = V- (Av) -0a(*v). 



(91) 



where e = 1,0 for the exact dynamics (ED) or ZA, re- 
spectively, we have used the fact that 



V- [(vV)v] = £ 2 ($ v )+V-(v# v ), 



(92) 



and that the ZA corresponds to replacing the Poisson 
equation by V 2 $ = 3H 2 n m 9/2 [71 [75]. 
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Equations ( 90||9~i" | make clear that the evolution of \ 
depends only on the velocity field through A. Further- 
more, by construction, x an d A vanish in linear theory, 
and since matter fluctuations are i n th e growing mode, 
A(") oc e"*, and one can solve Eq. i91l immediately 



A (n) 



2n + 3e 



[V ■ ^ 



(93) 



This gives a recursion relation starting from AW = 0. 
Note that denotes the n th -order contribution to Qi 
because $ v is a non-linear quantity. The second-order 
solution is straightforward, 



A (2) 



3e 



^ 2) ($v). 



(94) 



When used in Eq. (90), this gives (upon integration) the 
desired second-order solution: 



A 2 ) 



|[^] 2 - 2 f^(l-e^)^($ v ), (95) 



after using the initial condition x* 2 ' = (62/2) [&* 1 ''] 2 and 

(96) 



6 2 = b* 2 e~ 2y . 



Comparing Eqs. (89) and (95) we then identify 



72 



2(6! 



1) 



4 + 3e 



(1 



(97) 



Equation (95) agrees with Eq. (43) for the ED (e = 1), 
after recalling that the structure of Q2 yields monopole 
and quadrupole equal to each other but opposite in sign. 
We see from this that the change in dynamics (ED or 
ZA) only affects the amplitude of 72, not the structure 
of the second-order solution. 

Let us now find out the third-order solution. From 



Eq. ( 90 ) we have after using Eq. ( 93 ) 



X (3) = fc^[V.(A( 2 )v( 1 ))-^($ v )] + V.(x^v( 1 )), 

which can be integrated to give, after using the in itia l 
condition X i 3) = b* S^S^ + {b* 3 /6) $ 1} ] 3 , and Eqs. {94} 
and 1)5]) 



x < 3 > = M (1) <5 (2) + ! [<5 (1) ] 3 - 



2 6 2 



2(6! - 1) 
4 + 3e 



6 

8 + 3e 
6(2 + e) 



+ 



3e 
3e + 4 
6(2 + e)' 



V • 



(^ 2) ($v)v^ 



) - 



61-1 

6 + 3e 



(l-e-^)^ 3) ($ v ), 



(99) 



where we used 5^ = [S^]' 
and 



.(i).V^ 1 



-2^ 2) /(4- 



■3e) 



63 = 63 e 



-3y 



36 2 (l-e-«). 



(100) 



The first two terms in the first line of Eq. ( 99 ) are pre- 



cisely those expected of local bias; the last is the mix- 
ing term between local and non-local bias induced by a 



non-zero 63, which identifies (3 — 62/(61 — 1) in Eq. (89) 
independent of the dynamics. 



To extract 73 from Eq. ( 89 ) , and deviations from this 



equation in going from the ZA to the ED, we must first 
subtract from the second line of Eq. fl99"l ) the third-order 

contribution 72 implied in Eq. (I89J). This gives, for 
the only remaining non-local contribution to third order, 



rNloc 



61-I 



(4 + 3e)(6 + 3e) 



+ 3s) - 6(2 + e)e- y + (4 + 3e) e - 2 » ^ Jj ($ v ) - V • {Q^ 1 ($ v ) vW) 



?(3)/ 



7(2), 



(101) 



r 



To compute Q\ ; ($ v ), we need the non-linear evolution 
of the velocity potential to second order, $ v = +<I>i 2 ' , 
which gives 



ha 



52 2) (1>v) = (V. 13 $«) 2 -(V 2 $( 1 

^ 3) ($ v ) = 2(V ii H 1) Vi i $i 2) - V 2 $WV 2 $| 



(102) 



2h 



(103) 



Our calculation so far holds for both the ZA and the 
ED. While ' only depends on the linear potential (and 

thus it is independent of dynamics), Q\ depends on the 

(2) 

details of the dyna mics through <H ■ Therefore, we first 
compute Eq. ( 101 1 in the ZA, for which the second-order 

,(2) 
,ZA 



potential is straightforward: 



[V I $l 1) ] 2 /2, and 
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using Eqs. p2p03l ) we see that, remarkably, 



0$a(*v) - v • (^ 2) ($ v ) v«) = g 3 ($ v ). (104) 

This is the third Galileon operator to leading order. 
Therefore, this completes the proof of Eq. (89), with a 



73 coefficient that can be obtained from the amplitude in 
Eq. 1101b by setting e = 0. 



We now extend this result to the ED by using the 
fact that the ZA gives rise to the exact Galileon £3 in 



Eq. (104). We can obtain the ED second-order potential 



by noting that 6^ 2) = [£«] 2 +vW • ; /( 4 + 3e )> 
and thus 



y2 $ (2) = V 2 $ (2) 



ZA 



3e n {2) 



(4 + 3e) 



g?>(* v ). (105) 



From this we deduce 



3e 



$(2) _ <j>(2) 1 " $, IPi 

v ' ZA+ (4 + 3e) 2LPT ' 



(106) 



where the potential ( &2LPT is precisely the 2LPT poten- 
tial for the displacement field to second-order. It obeys 
the Poisson equation [49l [76] 



(107) 



V 2 $ 2LPT - -G 2 H% 



which implies that Eq. ( 104 ) becomes 



# ($v)-v-(0f (* v )vW) = g,(* v ; 

6e 

+ 



T(2)/ 



(108) 

(4 + 3£) 2 (^,$ 2L px) 
where, with a slight abuse of notation, we have defined 

OipP, $2LPT) = Vy$WVy$2LPT ~ V 2 $«V 2 $ 2L p T . 

(10?) 

We thus obtain the main result of this section, that in 
the exact dynamics the non-local part of the bias relation 
is given by (setting e — 1) 



rNlo 



I2Q2 (S v )(l + /3<f) 



6 



(110) 



+ 73 ( £ 3 ($v) + ^2($t 1) ,*2LPT; 



with coefficients given in terms of the local bias parame- 
ters by, 

72 = - 1)(1 " e""), - ^- &2 - I y, (111) 



73 = - 1)(1 - e -")(ll - 7e-»). 



(112) 



Equation (110) replaces Eq. (89) and includes the fact 



that the dynamics of gravitational instability is non-local; 
as a result of this, the bias relation is not only non-local 
in terms of the density, but also in terms of the velocity 
potential, through the non-local dependence o f $2J, pt on 



$v from inverting the Poisson equation Eq. (Il07|). This 



term can also be written using the (scaled) gravitational 
potential $, 



7 r 



Gi&P, $2LPT) = T £2 (*) " G2 (*v) 



(113) 



where V 2 <£> = 8. Therefore, one expects that the non- 
local bias relation to any order can be written in terms 
of Q2 ($ v )i Q2 ($), G3 ($v), and G 3 ($) including nonlin- 
ear combinations of them and multiplications of them by 
powers of 6. 

In the discussion so far we have assumed for simplicity 
that bias at formation is local. If we allow for non-local 



bias at formation, adding to Eq. (82) non-local terms 
given by 



e*Nloc 



72 &(«;) + 73 ea(*v) 



(114) 



it is easy to check that Eq. ( 110 1 still holds, but with the 
coefficients in Eqs. ( 111|112 ) changed to 



73 -> 73 



,-21/ 

72*e- 2y (l 



0+73' 



-3y 



(115) 
(116) 



In Appendix [A"| we generalize these results deriving how 
the 6™ oc obtained here can be related to that in the case 
when tracers are not conserved. 

Based on symmetry arguments 34 suggested a model 
of non-local bias with similar, but not identical, struc- 
ture to what is derived here. Our results in the absence 
of velocity bias, Eq. (110), agree with them to second- 
order, while at third order our results differ somewhat: 
their Ss 2 term corresponds to our term proportional to 
/3, their s 3 term can account partially for Q3 and their st 
term basically corresponds to our non-local term given 
by Eq. (113). However, they include an extra term to 



third order (which they call iff) which we find unneces- 
sary as it can be written in terms of the other third-order 
contributions. Thus we find that the most general third- 
order non-local contributions in the absence of velocity 
bias contain only three, not four, free parameters. In our 
case, these parameters are not free as they are connected 
by dynamical evolution to local bias (Eq s. |111||112 ), al- 
though if formation bias is non-local (Eq. 114) there are 
extra free parameters (Eqs. 115||116| ). 



VI. 



LARGE-SCALE NON-LOCAL BIAS IN 
NUMERICAL SIMULATIONS 



We now discuss how the ideas presented in the pre- 
vious sections can be used to ascertain the extent to 
which the large scale bias of dark matter halos in nu- 
merical simulations is local. We use a set of 50 simula- 
tions, each containing N par — 640 3 particles within a co- 
moving box-size of side Lbox = 1280 ft. -1 Mpc. The total 
comoving volume is thus approximately 105 (h~ 1 Gpc) 3 . 
Cosmological parameters were ft m — 0.27, 17a = 0.73, 
fib = 0.046 and h = 0.72, together with scalar spectral 



16 



TABLE I. Halo Samples used in this paper 



Halo Sample 


z 


bx 


Mass bin [10 13 M Q //i] 


LMzO 





1.43 


4-7 


MMzO 





1.75 


7-15 


HMzO 





2.66 


> 15 


LMz0.5 


0.5 


1.88 


3-5 


MMz0.5 


0.5 


2.26 


5-10 


HMz0.5 


0.5 


3.29 


> 10 


LMzl 


1 


2.43 


2 - 3.1 


MMzl 


1 


2.86 


3.1 - 5.7 


HMzl 


1 


3.99 


> 5.7 



index n s — 1 and normalization ag = 0.9. The simula- 
tions were run using Gadget2 [77] with initial conditions 
set at Zi = 49 using 2nd order-Lagrangian Perturbation 
Theory (2LPT) @9l [78]. The halos are identified using 
the friends-of-friends algorithm with linking length equal 
to 0.2 times the mean inter-particle separation. We di- 
vide our halo sample into three mass bins at each redshift 
z = 0,0.5,1. Table [I] shows the main features of each 
of these, including the large-scale (linear) bias obtained 
from measuring the cross-power spectrum between halos 
and matter, i.e. b x = Ph. m / Pmm, and averaging over 
scales k < 0.05 ft Mpc -1 . 

To assess the locality of large-scale halo bias in the sim- 
ulations we proceed as follows. We build the smoothed 
matter fluctuations by interpolating the dark matter par- 
ticles in the simulation to a grid of size -ZV gr id = 180 (cor- 
responding to a grid separation of ~ 7 h^ 1 Mpc), Fourier 
transforming using FFT's, multiplying by the Fourier 
transform of a real-space top-hat window function of ra- 
dius R s — 40 /i -1 Mpc, and Fourier transforming back to 
real space. We build the smoothed halo overdensity field 
similarly. We build the smoothed Galileon fields Q 2 and 
C?3 from the velocity field by first constructing the velocity 
potential $ v (and velocity divergence 9 V = V 2 $ v ) by us- 
ing a Delaunay tessellation to build the volume weighted 
velocity field on the grid (see [75] for details), construct- 
ing the Galileon fields on the grid and then smoothing 
them as one does for any scalar field (S or S^) as ex- 
plained above. That is, 

g 2 (x) = J e - ik ^ x (/x? 2 - i)9 v {k 1 )8 v (k 2 )w 12 d 3 k 1 d 3 k 2 

where W\ 2 = W{k\ 2 R s \ Hij = k\ ■ kj and 

Ga(x) = J e- lk ™ x (1 + 2/1x2^3^31 - M?2 - l4a " Mil) 
x 9 v (k 1 )9 v (k 2 )8 v (k 3 ) W 123 d 3 kid 3 k 2 d 3 k 3 . (118) 

We ignore, for simplicity, the extra non-local term de- 
pending on the 2LPT potential (see Eq. 
this term into the plots we present in t 
not change the results. 




110). Including 



ris section does 



FIG. 4. Illustration of non-local large-scale bias in numer- 
ical simulations for high-mass halos at z = 1 (see HMzl 
in Table El. The plot shows surfaces of constant Sh = 
—0.3,0.1,0.5,0.9 (from left to right, or red, blue, yellow, and 
green, respectively) as a function 5, Q2 and C/3- If large-scale 
bias were a local function of S, surfaces of constant Sh would be 
8 — const, planes (see next figure). Instead, there is significant 
tilt (VSh is not parallel to the 5-axis) showing a non-negligible 
dependence on <?2- All fields (5, Q2, Qz and 5h) have been 
smoothed with a top-hat window of radius R s — 40 hT 1 Mpc. 



Note that since the Galileon fields are non-linear com- 
binations of (derivatives of) the velocity potential, this 
procedure is not the same as building the Galileon fields 
of the smoothed velocity potential, which would remove 
mode- couplings of the smoothing scale to smaller scales. 
This means that our smoothed Galileon fields depend to 
some extent on the choice of grid size (which effectively 
determines up to what scale we allow mode-couplings; in 
our case this is down to ~ 7 h^ 1 Mpc). However, since the 
velocity power spectrum is suppressed compared to the 
density at small scales [80], the dependence is not very 
strong, particularly because, in Q 2 , the coupling to small- 
scale modes requires wave vectors to be anti-colinear in 
which case their contribution to Q 2 vanishes. We have 
studied what happens if we increase N gri d and we see no 
significant change to the results presented below except 
for an increase in noise (from coupling to even smaller- 
scale modes). This is somewhat expected as one starts 
to probe couplings to scales comparable or smaller than 
the Lagrangian size halos. Ideally, one would use a grid 
size different for each halo sample so only scales larger 
than the respective Lagrangian radius are included in 
Eqs. pTfpTS l. 



As a result of this procedure, at each grid point in the 
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FIG. 5. Same as Fig. [4] but for low-mass halos at z — (see 
LMzO in Table [l|. For the least biased objects in our samples, 
bias becomes local. 



simulation box we have four fields smoothed at large- 
scales (R s = 40 hr 1 Mpc): 5, Q 2 , Gz and Sh- If large- 
scale bias were local, Sh would depend only on <5, and 
thus surfaces of constant Sh should agree with those of 
constant S, in other words, VSh in this three-dimensional 
space (5, G2, G3) should be parallel to the 5-&xis. Figure|4] 
shows this construction for the highest mass bin at z = 1, 
where the effects of non-local bias are the strongest: there 
is a clear tilt of the surfaces of constant Sh in the G2 direc- 
tion, but no discernible dependence on Gs- Therefore, in 
cells of fixed S, where local bias would predict a constant 
Sh, we see significant variations in Sh that scale with the 
value of Qi ■ Note that at fixed S, Sh is a decreasing func- 
tion of G2, as predicted by our simple arguments in the 
previous sections. 

Figure [5] shows what happens in the lowest-mass bin at 
z = 0, for which the large-scale linear bias is the smallest 
among our halo samples. We see now that bias does be- 
come local: planes of constant S agree with constant Sh, 
and V<5/j points along the 5-axis. This is also in qualita- 
tive agreement with our simple model developed in the 
previous sections. 

In Figure [6] we show the same plot for high-mass halos 
at z = 0, demonstrating that indeed more biased objects 
at fixed z do show more non-local large-scale bias. Again, 
\7Sh has a significant component in the ^-direction, and 
little (if anything) in the C/3-direction, showing that the 
results presented in Fig. [4] are generic. Our model in 
the previous section does predict the dependence on G2 
to be stronger than that on G3 but only by a factor of 
about two or so (the precise value depends on "forma- 
tion" time) ; the results from the simulations suggest that 




0.4 

FIG. 6. Same as Fig. [4] but for high-mass halos at z = (see 
HMzO in Table 



the suppression of the G3 amplitude is even greater. 

It is rather common (see e.g. [42, 81 83 for recent ex- 
amples) to present the bias relation from simulations in 
terms of a scatter plot of Sh and S, which corresponds to 
projecting out the G2 and G3 directions in our Figs. [4]- 
[6] Because of the tilt in the Q2 direction, a bias that is 
completely deterministic in 5, Q2 and G3 will lead, when 
projected into the <5-axis, to a stochastic Sh vs 8 relation 
with the scatter simply coming from points with the same 
S that have different Q 2 . The question which arises is if 
the scatter seen in the Sh vs S relation can be explained 
by this projection effect, at least partially? One way 
to address this is to ask whether the scatter about the 
tilted planes with constant Sh in the three-dimensional 
space (5, G2, G3) is significantly less than that seen in the 
ID scatter plot of Sh vs S. We find that indeed the mul- 
tidimensional scatter is smaller than the ID scatter, but 
only marginally so (with one exception, which we discuss 
in the next paragraph). This indicates that most of the 
scatter of the Sh vs S relation is not due to the depen- 
dence of Sh on the "hidden variables" G2 and £3. In fact, 
this scatter can be explained [2H [55] in the context of 
the excursion-set model of halo formation by noting that 
the small-scale density field (whose excursions above the 
collapse threshold correspond to halo formation) has a 
stochastic relation to the large-scale density field S. 

Having seen that there is little, if any, non-locality 
coming from G3 we look for the possible effects of ve- 
locity bias. From our model we expect that if there is 
velocity bias at the smoothing scale we consider (R s = 
40 h^ 1 Mpc), then a dipole non-local term V will appear 
in the bias relation. As discussed before, a statistical 
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FIG. 7. Same as Fig. [4] (high-mass halos at z = 1) but using 
the dipole field T> rather than Qz as the third axis. There 
is significant dependence not only on Qi but also on P, a 
signature of velocity bias, as expected for the most significant 
biased objects. 
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FIG. 8. Mean and dispersion of C/2 (top) and C/3 (bottom) 
as a function of 5 at z — for top-hat smoothing with R a = 
40 h' 1 Mpc. Solid lines show the expected G2 = — (2/3)5 2 and 
G:i = 0. 



velocity bias of halos is expected for the most biased ob- 
jects [551 [SB], an d while this statistical effect vanishes at 
large scales as fc 2 , the scale at which this is negligible 
can be very large for very massive halos. We thus repeat 
the process discussed above to construct T> smoothed on 
R s = 40 h^ 1 Mpc at each grid point, 

V(x) = J e ~ ik ™- x /iu(^) 0v(fci)0v(fc2)Wi 2 d 3 k lC l 3 k 2 

(H9) 

and replace the C/3-axis by a 2?-axis for the halos in 
Fig. [4] which are the most biased objects in our sam- 
ple. These results are shown in Fig. [7j and confirm there 
is a significant component of in the Indirection. 
This dipole dependence, unlike that on CJ 2 , quickly dis- 
appears as we consider less extreme objects, so it is rele- 
vant only for the most rare halos. It is however interest- 
ing to note that for this sample of highly biased objects, 
the multidimensional scatter of Sh at fixed S,Q2,T>, e.g. 
5,02,1) = 0.15,0.08,0 is 0.35, substantially smaller than 
the ID scatter of 0.48 at fixed 5 — 0.15. This indicates 
that a significant fraction of the scatter in the Sh-S rela- 
tion may be due to projection of non-local bias for highly 
biased halos. 

Finally, it is useful to ask how the mean and dispersion 
of Q2 and C/3 vary as a function of 5. In Figure [8] we plot 
the mean (symbols) and dispersion (error bars) of G2 (top 
panel) and C/3 (bottom panel) at fixed S as a function of 
8 at z = 0. The solid lines give the naive expectation 
corresponding to averaging the angular dependence of 
the Q 2 and G3 kernels, i.e. CJ 2 = — (2/3)S 2 and CJ 3 = 0, 



which are quite reasonable approximations. Therefore, a 
non-local bias relation of the form 

S h = h 6 + -1 S 2 + 72 G 2 + 73 G 3 (120) 

will look, when projected into the <5-axis, as an effective 
local bias with 

sf^s+ Q*- 472/3) S 2 (m) 

The relation &| = b 2 — 472/3 for the effective local 
quadratic bias will be useful to interpret the difference 
in quadratic bias parameters between local vs non-local 
fits to the bispectrum obtained below. 

We now consider the effects of non-local bias from a 
different angle, in how it impacts the bispectrum. Along 
the way we quantify more precisely the amplitude of non- 
local bias terms at large scales. 

VII. NON-LOCAL BIAS AND THE 
BISPECTRUM 

The bispectrum is sensitive to the actual couplings at 
second-order in perturbations in both the dark matter 
and the bias relation. Therefore, it can give a precise 
determination of how non-local bias is, in particular the 
amplitude of the CJ 2 dependencies found from smoothed 
fields in the previous section. Here we concentrate on 
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the bias relation up to second order, we thus include only 
local quadratic bias b 2 and the amplitude of the non-local 
effect through Q 2 . As we found in the previous section, 
there is no significant detection of a Q 3 dependence, and 
the (quadratic) dipole dependence is only significant for 
the most biased samples, which correspond to extremely 
rare halos. Thus we study the bias relation given by 



Sh = bi S + — S 2 + 72 Q 



(122) 



where we recall that for 72 = — 2(&i — l)/7, this corre- 
sponds to assuming local Lagrangian bias (ie. y — > 00 
in our simple model). In this way we can simultane- 
ously test for local Eulerian (7 = 0), local Lagrangian 
(7 = — 2(6i — l)/7) and more generic non-local bias. We 
note that a similar test (of local Eulerian vs Lagrangian 
bias) was performed in the PSCz galaxy survey bispec- 
trum [3] with the result that Eulerian local bias was a 
slightly better fit to the galaxy bispectrum. Our tests in 
this paper are in a very different regime, as PSCz galax- 
ies are anti-biased while our halos are positively biased 
(see Table |l|. 

To avoid dealing with the complications of the inac- 
curacy of Poisson shot-noise subtraction for halos (see 
appendix A in [ST]) that complicates interpreting devia- 
tions from the local bias description, here we just study 
the halo-matter-matter bispectrum bh mm rather than the 
the halo bispectrum [531 HI] > 



(S h (k 1 )S(k 2 )S(k 3 )) = S B {k 123 ) b hmm {k 1 ,k 2 , k 3 ). (123) 

Note that this is not a symmetric function of the wavec- 
tors. We thus define a symmetrized quantity [83 , 



-Bhmm — (bhmm ~t~ ^mhm ± ^mmh) j (-^^) 



which from Eq. ( 122 ) obeys 



B 



hmm 



= h B l 



23 



Si 



23 



123 



(125) 



where -B123 is the matter bispectrum and 

£123 = P1P2 + eye, /C123 = (M12 - 1) Pi P2 + eye, 

(126) 

with \i\ 2 the cosine of the angle between k\ and k 2 . The 
kernel K,\ 23 vanishes for colinear trianges where fiij = 
±1, thus the non-local correction is most important for 
isosceles triangles. 

We measured the halo-matter-matter bispectrum 
Bhmm and matter bispectrum Bi 23 from the 50 realiza- 
tions of the simulations at the three redshift outputs. 
The triangles included in the bispectrum analysis corre- 
spond to all triangles with sides from twice the fundamen- 
tal mode (2k f ~ 0.01 /iMpc -1 ) up to k < O.l/iMpc -1 , 



TABLE II. Local Eulerian bias parameters b\ and 62 obtained 
from halo-matter-matter bispectrum fits for all triangles with 
k < 0.1 /iMpc -1 . We also include the large-scale bias b x ob- 
tained from the halo-matter power spectrum, to be compared 
with 61. The last column indicates the goodness of the fit 
assuming a diagonal covariance matrix (iVdof = 148). 



Smrmlp 


ux 




bn 


v 2 /dof 


LMzO 


1.43 


1.42 ±0.01 


-0.91 ±0.03 


1.86 


MMzO 


1.75 


1.71 ± 0.01 


-0.55 ±0.03 


1.29 


HMzO 


2.66 


2.37 ± 0.02 


2.98 ±0.07 


3.74 


LMz0.5 


1.88 


1.77 ±0.01 


-0.15 ±0.03 


0.91 


MMz0.5 


2.26 


2.13 ±0.01 


0.67 ±0.03 


0.87 


HMz0.5 


3.29 


2.84 ± 0.03 


5.89 ±0.10 


3.77 


LMzl 


2.43 


2.22 ±0.01 


1.27 ±0.04 


0.89 


MMzl 


2.86 


2.62 ±0.02 


2.77 ±0.06 


1.07 


HMzl 


3.99 


3.41 ±0.05 


9.98 ± 0.14 


3.42 



binned in units of 2k f, yielding 150 binned triangles (cor- 
responding to ~ 1.2 x 10 8 fundamental triangles of sides 
inside the prescribed bins and all possible orientations). 
This together with the measured matter power spectrum 
can be used in Eq. (1251 to fit for the parameters 61, b 2 
and 72. In what follows we discuss such constraints for 
all the halo samples. 



Table [TT] shows the results from fitting Eulerian local 
bias (72 = 0) to the relation in Eq. (1251. For compari- 
son, in this and other tables, we reproduce the value of 
the large-scale linear bias obtained from the halo-matter 
cross spectrum b x . Note that for the lowest biased ob- 
jects in our sample, LMzO, the linear bias obtained from 
the bispectrum 61 agrees with b x , but this agreement 
disappears for all other samples, giving a significantly 
smaller bi than the large-scale linear bias b x shown by 
the power spectrum, increasingly so for more biased ob- 
jects. Recently, [53] found a similar result for halos with 
more than 20 particles at z = 0. Here, we highlight the 
mass and redshift dependence of this issue in more detail. 
As shown in [83] . had we used the reduced bispectrum 
Q = B/E rather the bispectrum itself to find the bias pa- 
rameters, then we would have found the opposite result, 
i.e. a linear bias 61 smaller than b x . To explain why, let 
us for definiteness define a reduced halo-matter-matter 
bispectrum by 



Q 



Bhr, 



B, 



hmm 



hmm 



(Px(fci)Px(fc 2 )±cyc.) & 2 X E 123 



(127) 



Thus, while the halo-matter-matter bispectrum fits yield 
b\B and b 2 B with Bhmm = b\sB ± b 2 B E/3, the reduced 
bispectrum yield parameters b\Q and b 2 Q with Qhmm = 
Q/biQ + b 2 Q/3b 2 Q. These are related by 



he 
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(128) 
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TABLE III. Eulerian bias parameters 61 and 02 obtained from 
doing a Lagrangian local bias model fit to the bispectrum. 



o£impl6 
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IjIVIZU 


1 A*\ 
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1 /is -un m 


i 9^ -t- n fM 


9 1 9 


IVIIVIZU 


1.(0 


1.01 m u.ui 


— 1.10 It U.UO 


1 .OD 


HMzO 


2.66 


2.59 ±0.02 


1.78 ±0.07 


2.73 


LMz0.5 


1.88 


1.87 ±0.01 


-0.79 ±0.04 


0.94 


MMz0.5 


2.26 


2.30 ±0.01 


-0.26 ± 0.04 


0.72 


HMz0.5 


3.29 


3.12 ±0.03 


4.34 ± 0.11 


2.91 


LMzf 


2.43 


2.40 ±0.02 


0.27 ±0.05 


0.77 


MMzl 


2.86 


2.85 ±0.02 


1.45 ±0.06 


0.82 


HMzl 


3.99 


3.77 ±0.05 


7.97 ±0.16 


2.74 



TABLE IV. Eulerian bias parameters &i and 62 and non-local 
72 parameter obtained from doing a quadratic non-local bias 
model fit to the bispectrum. For comparison purposes, note 
that a non-zero 72 gives an effective —(4/3)72 contribution to 
62 (see top panel in Fig. |8j. Here N doi = 147. 



therefore, if bis > b x , then b\Q < b x . Similarly, for 
halo bispectra (rather than halo-matter-matter), the re- 
lationship between reduced and im-reduced bispectra lin- 
ear bias is instead &iq = b x (fe x /&is) 3 , an even bigger 
difference (i.e. the relative deviation of biq from b x is 
three times larger than for bis)- These disagreements 
will be resolved shortly by including non-local bias. 

Table |III| shows the analogous results when the bias 
is assumed to be local in Lagrangian space, equivalent 
to assuming 72 = -2(b x - l)/7 in Eq. ( [X22"] ). The re- 
sults in this case are somewhat mixed. At z — the 
results are worse than for the Eulerian case, except at 
high mass. At higher redshifts, the Lagrangian results 
show improvement, particularly at z = 1, but there are 
still some significant discrepancies between b\ and b x , 
and in any case the x 2 /dof are not very convincing. 

Finally, Table |IV| shows the results for the non-local 
bias model with the amplitude of Qi being fit for. The 
results show now a significant improvement, in particular 
b\ is always within two-sigma of the b x values, for all 
redshifts and halo masses considered. We note that the 
average (over all halo samples) deviations of b\ from b x 
are He, 4.5<r and 1.5a for Eulerian, Lagrangian and non- 
local bias fits to the bispectrum, respectively. Thus we 
reject local Eulerian and Lagrangian bias models at high 
significance. The price to pay in fitting for 72 as well is 
an increase in the b\ error bars, by a factor of almost two. 

The values for 72 in Table [TV| show a clear dependence 
with linear bias, which is plotted in Fig. [9] (using the 
more precise value of b x as linear bias). We see that the 
results fall mostly along along a "universal" line given by 
— 2(6! — 1.43)/7 (solid line), except for the most biased ha- 
los at each redshift which fall below this line (closer to the 
Lagrangian bias result, shown in dashed line). However 
it is precisely these highly biased objects that may have 
extra non-local contributions (such as a dipole, as dis- 
cussed in the last section), so it is not clear at this point 
how reliable this behavior is. On the other hand, note 
that the solid line in Fig. [9] is not a fit to the data, but it 
serves to illustrate deviations from local Lagrangian bias 
for our least biased samples. More work is needed to see 



Sample 


Ox 


61 


02 


72 


X 2 /dof 


LMzO 


1.43 


1.42 ± 0.02 


-0.92 ± 0.08 


-0.01 ± 0.03 


1.87 


MMzO 


1.75 


1.76 ±0.02 


-0.81 ±0.08 


-0.10 ±0.03 


1.19 


HMzO 


2.66 


2.61 ± 0.04 


1.71 ± 0.18 


-0.48 ±0.06 


2.74 


LMz0.5 


1.88 


1.83 ±0.02 


-0.46 ±0.09 


-0.12 ±0.03 


0.84 


MMz0.5 


2.26 


2.24 ± 0.02 


0.05 ±0.09 


-0.24 ±0.03 


0.67 


HMz0.5 


3.29 


3.16 ±0.06 


4.10 ± 0.28 


-0.70 ±0.10 


2.91 


LMzl 


2.43 


2.35 ± 0.03 


0.57 ±0.13 


-0.28 ±0.05 


0.74 


MMzl 


2.86 


2.80 ±0.03 


1.70 ±0.16 


-0.42 ±0.06 


0.80 


HMzl 


3.99 


3.84 ±0.08 


7.55 ±0.41 


-0.96 ±0.16 


2.73 



whether one could understand these results from theoret- 
ical arguments. We note however that it is not surprising 
that bias is not local in Lagrangian space, even in simple 
extensions of the excursion set of halo formation the bar- 
rier for collapse is known to depend on other quantities 
than the overdensity <5, mostly on the ellipticity param- 
eter e [32]. Appendix [B] shows the relationship between 
ellipticity e, prolateness p and the invariants of the de- 
formation tensor or Galileons. 

We see then that the presence of non-local bias (G2) re- 
quired from the multi-dimensional plots in the previous 
section is confirmed by the bispectrum analysis, which 
shows that including such terms solves a systematic er- 
ror in the determination of the linear bias, increasing 
for more biased objects. This is important because this 
systematic error would otherwise affect the determina- 
tion of cosmological parameters from a bispectrum anal- 
ysis (see |85] ) , particularly for luminous galaxies (such as 
LRGs in SDSS) that populate high-mass halos. The extra 
dependence on Q2 is also important in at least two more 
aspects: it introduces a dependence on triangle shape 
that is degenerate with brane-induced modifications of 
gravity [SHI HZ] , and also mimics an equilateral- type pri- 
mordial non-Gaussianity signature (see Fig. 1 in 88 ). 
Therefore, for all these reasons, it is important that such 
dependencies are taken into account when doing bispec- 
trum analyses in galaxy surveys, extending what was 
done already in [5] by considering both Eulerian and La- 
grangian local bias models. 



VIII. CONCLUSIONS 

In this paper we studied the non-localities induced 
in the bias relation by gravitational evolution, provid- 
ing results under a number of different scenarios. In the 
simplest case, galaxies form at a single time and evolve 
conserving their comoving number density (no merging) 
following the dark matter (no velocity bias). In this 
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0.2 




FIG. 9. Dependence of the Q2 non-local amplitude 72 (see 
Table |TV| ) bispectrum fits on the linear bias (61 = 6 X ) for the 
three halo mass bins at three different redshifts, z = (blue 
triangles), z = 0.5 (red squares) and z = 1 (green pentagons). 
The dashed line shows the predictions of local Lagrangian 
bias, and the solid line corresponds to 72 = — 2(fci — 1.43)/7. 



case, even if galaxies initially are locally and linearly 
biased, they develop non-local and non-linear bias that 
can be best described by the invariants of the deforma- 
tion tensor (Galileons, see Eqs. 84p5 ). The main signa- 
ture of this is a quadrupole field to second-order (already 
present in [35j, and best known from local Lagrangian 
bias [301 SI])- We also derive for the first time the evolu- 
tion of the bias to third order (Eq. 110 1 and show that it 
is not only non-local in the density but also in the poten- 
tial, as a result of the non-local nature of gravitational 
instability. 

We generalized the second-order result in several direc- 
tions. First, we considered what happens when merging 
and formation of galaxies is taken into account, show- 
ing that the non-local bias to second-order can still be 
written due to symmetry reasons in the same Galileon 
form but with an amplitude that depends on the evo- 
lution of the comoving number density that cannot be 
made to vanish (and thus keep bias local) unless galaxies 
are unbiased. Along the way, we generalized the stan- 
dard formula of the time evolution of linear bias I2S1 EH 



to the case where merging is important (see Eq. 73 1 . We 



also extended the results to the case where there is ve- 
locity bias b vi finding its time evolution in linear theory, 
Eq. ( 32 ) . We showed that when there is velocity bias then 



an extra non-local contribution appears at second-order. 
This is related to the breaking of Galilean invariance in 
the bias relation proportional to (b v — 1). This relative 
motion between tracers and matter generates a dipole 



term (see Eqs. 37p9 ). Appendix |A| generalizes these re- 
sults to the case when tracers are not conserved, finding 
how the non-local bias relation in this case is related to 
conserved tracers, even if bias at formation and merging 
is non-local. 

To test these results we proceeded in two different 
ways: by studying the bias relation in simulations at the 
leve l of the halo field compared to the matter field (Sec- 
tion VI I , and statistically by measuring the halo bispec- 



trum in simulations and comparing to the predictions of 
non-local bias (Section |vTi|). 



Regarding the halo bias relation, we found that halo 
ovcrdensities in cells of fixed dark matter density indeed 
vary with the value of the quadratic Galileon field (de- 
scribing the strength of the tidal field) as expected from 
our calculation (with stronger dependence for more bi- 
ased objects, see Figs. [4 [6]). This is a direct demonstra- 
tion that such non-loca terms are present in the bias 
relation for dark matter halos. For highly biased halos 
at high redshift we also found some evidence for the ef- 
fects of a dipole field (see Fig. [7]), which might be related 
to velocity bias. More work is needed to establish this 
precisely in the simulations. 

The halo bispectrum measurements in the simulations 
show a clear detection of the effects of non-local bias as 
the linear bias increases (see Table IV and Fig. [9]) , fix- 



ing a systematic error in the determination of linear bias 
(also seen recently in |83| ) when assuming local Eulerian 
bias (see Table [n ) . We also performed local Lagrangian 
bias fits to the ha o bispectrum, showing that this model 
also leads to systematic errors in the determination of 
the linear bias, although it becomes better for highly bi- 
ased objects (but the resulting chi-square values are still 
large). We obtain an approximate relation valid for the 
least biased halos in our sample (solid line in Fig. [9 ) be- 
tween the amplitude of the non-local term 72 and the 
linear bias inspired by a simple modification of the local 
Lagrangian bias model that would be interesting to test 
in other cosmologies and in other linear bias regimes. Our 
results indicate that viewed in Lagrangian space bias is 
also non-local, this can be checked simply by redoing the 
multidimensional scatter plots presented in Section [VT| in 
Lagrangian space. 

Accounting for these effects when modeling galaxy bias 
is important for correctly describing the dependence on 
triangle shape of the galaxy bispectrum, and hence con- 
straining cosmological parameters and primordial non- 
Gaussianity. The non-locality of bias also introduces new 
angular dependences in the redshift-space power spec- 
trum and bispectrum (since it modifies the angular de- 
pendence of perturbation theory kernels) that will be im- 
portant to include and test against simulations. We hope 
to report on this in the near future. 
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Appendix A: More on Non-Conserved Tracers 

We have demonstrated in Section IIVI that even when 
the tracers are not conserved, due to formation and merg- 
ing, the structure of the non-local bias up to second or- 
der is the same as in the conserved case. Here we extend 
these results to third order and beyond. Furthermore, 
here we also keep the source function j general, allowing 
a dependence on other fields in addition to 8, allowing for 
non- locality in the galaxy formation/merging processes. 

We write rig = Jig (1 + S g ) and introduce a generic 
dependence of the source on fluctuations through ip by 



j( P ,e)=j(p,Q) [i+mo] 



(Al) 



where, in addition to (5, the source function j also de- 
pends on 0, which collectively denotes any extra fields 
(may even be non-local, e.g. Q2). More specifically, the 
fluctuating part of j is given by 



where 
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Note that if = (e.g. as in the case = Q2) one 
can still apply this, as c*( l is well-behaved in 0. Using 
Eq. JAl) in Eq. Q, we get 



d_ 

Oy 



(n g c) S g ) - 4 C) V • [(1 + S g ) v] = fiW i/,(5, C), (A4) 



where a dot denotes a derivative with respect to y and 
we have used Eq. (651 to eliminate the background con- 



tributions. We now construct the equation of motion for 
X = S g — biS (see Eq. 28) for the non-conserved case. It 
obeys, 

d - (4 c) x) = [(&i- 1) A 



dy 



'^cr Aeons 



f V-(xv)]nW+ftW(^-6J5) 
h^^-blS), (A5) 



where Xcons denotes the solution when the tracers are 



conserved (see Eq. 90 ). In deriving Eq. ( A5 ), we have also 
used Eq. (A4), Eq. (74) and the continuity equation for 

- (c) 

dark matter. When tracers are conserved, n g — and 
thus Eq. ( A5 ) reduces to Eq. ( 90 ) . The general solution 



for x can then be written in terms of that in the conserved 
case, namely 



4 c) x = 



dy* n* Xc 



dn» (V - &i<0, (A6) 



which says that the solution for non-conserved tracers 
is essentially that for conserved tracers weighted by the 
evolution of the number density (first term) plus a term 
that depend on the sources that describe galaxy forma- 
tion and merging (ip) . If these sources are local functions 
of 5, then this extra term does not lead to any non-local 
contributions, and the non-locality of bias is precisely of 
the same form as in the conserved case (with slightly dif- 
ferent coefficients that depend on the evolution of the 
comoving number density). If, on the other hand, galaxy 
formation or merging depends on non-local functions of 8 
(e.g. the velocity divergence, the tidal field) then an ex- 
tra non-local contribution to galaxy bias gets generated 
by the second term in Eq. ( A6 1 . 



All this implies that the non-local part of the galaxy 
bias in the general case is related to that for conserved 
tracers by, 



j(c)^Nloc _ 



dy* n* cL cons 



dn* ip 
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or, integrating by parts, 



rNlo 



= 8, 



Nloc 
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(A8) 



This gives our most general expression for the non- 
local part of galaxy bias when formation and/or merging 
cannot be neglected. Although we have implicitly as- 
sumed there is no velocity bias, it's easy to check that 



Eq. (A8) is also valid when there is velocity bias, which 
only changes <5^J,° c ns through dipole terms but not the 
relationship given by Eq. ( A8 ) itself. 



Appendix B: relation of Q2 and Q3 to e and p 

The terms induced by gravity are present due to tidal 
fields, described by derivatives of the gravitational poten- 
tial and velocity potential (proportional to each other in 
linear theory) . To describe the shape of the gravitational 
potential it is common to introduce the ellipticity e and 
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prolateness p, defined from the eigenvalues A; (i — 1, 2, 3) 



of V i3 -$: 



wondered if the rotationally invariariant quantities, 

7i=Tr(D)=^Ai = 5, (B4) 



/a — A1A2 + A1A3 + A2A3, 
I 3 = Defc(D) = A, 



(B5) 
(B6) 



5 = Ai + A 2 + A3, 
Ai — A3 

Ai + A3 — 2A2 
P = 26 ' 



(Bl) 
(B2) 

(B3) 



This set of parameters is used in triaxial evolution mod- 
els of nonlinear structure formation [3TJ However, 
because e and p are ratios of the eigenvalues, it is not 
obvious that they are the best choice of parameters in 
a perturbative analysis. In particular, one might have 



are more relevant. When expressed in terms of (5, e,p) 
these are 



/ 2 = y[l-(3e 2 +p 2 )], 



/ 3 = |*(l~2p)[(l+p) 2 -9e 2 
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(B7) 

Since the Ij do not depend on taking ratios of the eigen- 
values, they, or other quantities built from them, have 
considerable appeal. One such combination is 



I 2 - 3/ 2 



2/f - 9/1/2 + 27/3 



Despite the appearance of I± in their definition, r and u 
are actually independent of I\. Moreover, they are pre- 
cisely the quantities which arise in a perturbative anal- 
ysis of the ellipsoidal collapse model: J\ and J2 of [89 
are our r 2 and 9u 3 respectively. Notice that e, p, r 2 and 
it 3 all vanish for a spherically symmetric perturbation 
(A 1= A 2 = A3=<V3). 
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